
figure;
%% Glucose plot of KM vs. absolute temperature
KMmat = monod_KM_glucose;

T = KMmat(:,1);
KM = KMmat(:,2)/1000;
KMerr = KMmat(:,3)/1000;

errorbar(1./(T+273.15), log(KM), KMerr./KM, 'bo-', 'linewidth', 1);
hold on;
%% Acetate plot of KM vs. absolute temperature
KMmat = monod_KM_acetate;

T = KMmat(:,1);
KM = KMmat(:,2);
KMerr = KMmat(:,3);

errorbar(1./(T+273.15), log(KM), KMerr./KM, 'ko-', 'linewidth', 1);

%% Succinate plot of KM vs. absolute temperature
KMmat = monod_KM_succinate;

T = KMmat(:,1);
KM = KMmat(:,2);
KMerr = KMmat(:,3);

errorbar(1./(T+273.15), log(KM), KMerr./KM, 'mo-', 'linewidth', 1);


%% LB plot of KM vs. absolute temperature
KMmat = monod_KM_LB;

T = KMmat(:,1);
KM = KMmat(:,2);
KMerr = KMmat(:,3);

errorbar(1./(T+273.15), log(KM), KMerr./KM, 'ro-', 'linewidth', 1);

%% Succinate plot of KM vs. absolute temperature
KMmat = monod_KM_NZCasePlus;

T = KMmat(:,1);
KM = KMmat(:,2);
KMerr = KMmat(:,3);

errorbar(1./(T+273.15), log(KM), KMerr./KM, 'go-', 'linewidth', 1);

set(gcf, 'Position', [0 0 400 300]);
set(gca, 'fontsize', 20);
xlabel('1/T (1/K)');
ylabel('log({\it K}_M)')
xlim([3.2E-3 3.37E-3])
box off

